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We study quantum spin systems with a proper combination of geometric frustration, spin-orbit 
coupling and ferromagnetism. We argue that such a system is likely to be in a chiral spin state, a 
fractional quantum Hall (FQH) state for bosonic spin degrees of freedom. The energy scale of the 
bosonic FQH state is of the same order as the spin-orbit coupling and ferromagnetism — overall 
much higher than the energy scale of FQH states in semiconductors. 



I. INTRODUCTION 

Landau symmetry breaking^ ^ has been the standard 
theoretical concept in the classification of phases and 
transitions between them. However, this theory turned 
out insufficient when the fractional quantum Hall (FQH) 
stat^^ was discovered. These states (FQH states and 
spin liquids) are not distinguished by their symmetries; 
instead they have new topological quantum numbers such 
as robust ground state degeneracjEHl and robust non- 
Abelian Berry's phases^. The topological ordei^ ^ asso- 
ciated with topological quantum numbers has been pro- 
posed for the classification of these states. Recently, it 
was realized that topological order can be interpreted as 
patterns of long range quantum entanglement^^ This 
long range entanglement has important applications for 
topological quantum computation: the robust ground 
state degeneracy can be used as quantum memorjEl 
fractional defects from the entangled states which carry 
fractional charges^ and fractional statistic^^JHMl ^q^. non- 
Abelian statistics^^ ^^) can perform fault tolerant quan- 
tum computation^^ 

Although it has attractive concepts and applications, 
topological order is only realized at very low tempera- 
tures in FQH systems^ ^. In this paper we present a 
proposal to realize highly entangled topological states at 
higher temperatures. The ideal is to combine geomet- 
ric frustration, spin-orbit coupling and ferromagnetism 
in quantum spin systems. Both spin-orbit coupling and 
ferromagnetism can have high energy scales and appear 
at room temperature. Their combination breaks time- 
reversal symmetry which leads to rich and complicated 
interference from quantum spin fluctuations. In this pa- 
per we show that they can lead to highly entangled topo- 
logical states at high temperatures. 

Quantum spins on the kagome lattice are ge- 
ometrically frustrated systems. They appear 
in the following compounds: Herbert smithite 
Zn Cu3 (OH)g Cl2,™t^ Kapellasite Cu3Zn(OH)gCl2,^ 
Yo.5Cao.5BaCo407F Mq,Cu^_,{0H)qC12^ 
CaBaCo^OjW^ Pr^Gar^SiOi^^ Nd^Ga^SiOi^W^ 
BaCus V2 Os (OH) 2 ,1^ Cu( 1 ,3-benzenedicarboxylate) ^ 
KFe3{0H)e{S0^^ YBaCo^Oj,^^^^ YBaCo^AlOj, 
YBaCo^FeOjW^ 7-^^12(0^)3^/,^^ Nir,(TeO?,)ABr2. 
7Vz5(Te03)4C/2,^ Cu3V20,{OH)2/2H20W^ 



Cs2Cu3CeFi2^ Cs2CusSnFi2, Rb2CusSnFi2F^ 
Cu2{OD)sClpi Cs2CusZrFi2, Cs2CusHfFi2^ and 
C03V2O8. Motivated by these materials, in this paper 
we study the Heisenberg model on the kagome lattice 
with additional spin-orbit interaction and Zeeman 
coupling ^^BzS^. Some related theoretical work 
can be found in Ref. I45|46l In Ref. |45] a model with 
spin-orbit interaction but no Zeeman coupling is studied; 
some mean-field spin liquid states are found. In Ref. 
46 a model with Zeeman coupling but no spin-orbit 
interaction is studied via numerical calculations. Two 
magnetization steps are found at M/Mmax = 1/3 
(stronger) and 2/3 (weaker) for a 36 spin cluster. 

In this paper, we study the state with magnetization 
(Sf) = 1/3. In section [nj we write down the quantum 
spin model with spin-orbit coupling on the Kagome lat- 
tice. In section [Hlj we map the spin model to the hard- 
core bosonic model in [III A| , construct three trial wave- 
functions for the polarized spin system (S'f) = 1/3 in 



IIIB and then evaluate the energy expectation for these 



three states in |IIIC[ We find that the bosonic quantum 
Hall state has the lowest energy. Lastly, we discuss the 



materials realization in HID In the Appendix [A[ we also 
discuss spin-orbit coupling in the transition metal oxide 
materials. 



II. QUANTUM SPINS ON THE KAGOME 
LATTICE WITH SPIN-ORBIT COUPLING 

The kagome lattice has 3 sites (labelled 1, 2 and 3) 
within every unit cell with the primitive vectors ai = 2ax 
and 3.2 = <^(x + V^y) {cL is the lattice constant), see 
Fig. [ija). The unit cell contains one hexagon and two 
triangles so it is geometrically frustrated. 

As shown in Fig. [l] (b), the triangle A123 on the 
kagome plane in Herbertsmithite ZnCu3(OH)gClJ2H^ 
contains three copper cations surrounded by distorted 
octahedrons sharing one chlorine corner while each pair 
shares an oxygen corner. Mediated by this oxygen, the 
Cu 3d^ electron hops from site ri to r2, e.g. see Fig. 

Inversion symmetry for Herbertsmithite breaks down 
explicitly, leading to a non-uniform charge distribution in 
the kagome lattice. For convenience, we model the charge 
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where D^j = wijD. Then 
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FIG. 1: (color onhne) (a) The kagome lattice with three dif- 
ferent sites / = 1, 2, 3 within the unit cell. Inversion symmetry 
breaking via a charge center in the hexagon leads to the ef- 
fective electric field j on the bond Yij , represented by green 
arrows which point from the middle of the bond to the cen- 
ter of every triangle on the kagome lattice, (b) The triangle 
Ai23 in the kagome lattice for Herbertsmithite^^"^. (c) The 
nearest neighbor bond Y12'. the electron hops from site ri to 
Y2 mediated by the oxygen atom. 



(4) 



Using standard second-order perturbation theory, we ob- 
tain the exchange term 



(5) 



Here J = At^ /U is the exchange coupling for the rotated 
spin operator = Qa^aa'C^^/. 

On the kagome lattice, we cannot find a gauge trans- 
formation as in Eq. ([3| that would be compatible for 
each site. So we have to write the Hamiltonian in terms 
of the original spin operators. On every bond, the ro- 
tated spin operators are related to the original ones as 
follows: 

Si = {I- cos{D)){hij ■ Si)hij + cos{D)Si 

— s'm{D)Si X hij (6) 

Sj = (1 — cos{D)){hij • Sj)hij + cos{D)Sj 

-^sm{D)Sj X hij (7) 

Thus we obtain the quantum spin model on the kagome 
lattice including spin-orbit coupling^^ 

H = J^{cos{2D)Si • S^- + sm{2D){Si x S^) • hij 

{ij) 

+2sin2(I))(S, •n,,)(5, -n,,)). (8) 



center as being in the hexagon, see FigjTJa). When hop- 
ping from ri to r2, the electron sees the electric field E12 
(labelled by the green arrow in Fig. [ija)). The effective 
electric field couples to the electron through the spin- 
orbit coupling vector D12 = Q^Ei2 x ri2 in the Rashba 
manner, where ri2 = ri — 

ti2 = X] (^exp(-i(J • T>i2)cjcj'c\^C2cj' + h.c.^ (1) 

(T(t' 

Here a = {(Jxi^yi^z) are the Pauli matrices. The coeffi- 
cient a should be chosen to make the spin-orbit coupling 
vector D12 dimensionless. Note that D12 = — D12. 

Including on-site interactions we obtain the Hubbard 
model with spin-orbit coupling for = 1/2 electrons on 
the kagome lattice 

H=-tY, ((e-*^-°-)<..'4c,>' +h.c.) + C/^nitn,;(2) 

ctct' i 

where i and j denote nearest neighbors. 

For a specified bond r^j, we can make a gauge 
transformation^^ 

c- ^c- = Vfe^^^/^^^-"^0 /c- / 

a' 

Cj(j y Cj(j ^ ^ (^ ^ ^ ^ ^ Cj(jf (3) 



III. POLARIZED SPIN STATE WITH 
TOPOLOGICAL ORDER 

A. Hardcore bosonic model 

For simplicity we choose the spin-orbit coupling vectors 
T>ij perpendicular to the kagome plane: D12 = D23 = 
D31 = Di/2' = D2/3/ = D3/1/ = Dz (only in-plane ef- 
fective electric fields E^j are considered). We use the 
Holstein-Primakoff transformation: 

s+ = bl ,s- = h ,S! = l-blbi (9) 

where bi is the hardcore bosonic operator = ^ij, 

= blbi < 1. This maps the spin model ^ onto a 
hardcore bosonic model^^ 

H = ^"^(^exp[{hij-z)i2D]blbj^h.c.^ 

(ij) 

-^J^riiTij (10) 

(ij) 

which describes interacting hardcore systems with hop- 
ping under effective fluxes as shown in Fig. [2] within the 
triangles A123 and Vi/2'3', there are fluxes = 6D; in 
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FIG. 2: (color online) Flux distribution for the hardcore 
bosons. 



the hexagon, there is flux (j)2 = — 20i. When 7^ 0,7r 
mod 27r, these effective fluxes break time-reversal sym- 
metry for this model. 

Now let us consider just one boson described by only 
the hopping term in the above Hamiltonian. The hopping 
Hamiltonian has three bands. We calculate the Berry 
curvatures over the Brillouin zone for the lowest band in 
the presence of the flux for different spin-orbit couplings: 
D = 0.025, 0.1 and 7r/8 . The Berry curvature is defined 
as follows 



Fn{k)=e,jdk^Af{k), AY"{k)=i{u„k\d,Junk) (11) 

where ii^k is the Bloch wave packet in the n-th band of 
the hopping Hamiltonian 



Ht\Un\^) = Cnkl^nk) (12) 



In this paper, instead of using kx and ky as the axes 
in /c-space, we use ki = kx and ^2 = (^x + V^ky)/2 
for convenience. The dispersion e^k has three bands (la- 
belled n = 6 for the bottom band, n = m for the middle 
band and n = t for the top band) as shown in Fig. [s] 
(a), (c) and (e). In all cases {D = 0.025, D = 0.1 and 
D = 7t/8) the three bands have nonzero Chern numbers 
Cb = 1, Cm = and Ct = —1, where the Chern number 

We plot the Berry curvature of the bottom band for 
D = 0.025, 0.1 and 7r/8 in Fig. |3](b), (d) and (f). We see 
that when D = 0.1, the lowest band is separated from 
the other bands by an energy gap and the lowest band is 
quite flat. Since the lowest band has a non-zero Chern 
number C5 = 1, it simulates the first Landau level in free 
space. By analogy to the quantum Hall effect in high 
magnetic field, the hardcore bosons are likely to form a 
u = 1/2 bosonic quantum Hall state when there is half a 
boson per unit cell. The boson filling number is / = 1/6 
per site which corresponds to the spin polarization 



{5f ) = 1/2 - / = 1/3 



(13) 



In other words, the polarized spin state (6'f) = 1/3 is 
likely to be a chiral spin liquid^ — a topologically or- 
dered state. 



(a). Band dispersions (D=0.025) (b). Bottom band curvature (D=0.025) 
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FIG. 3: (color online) (a), (c) and (e): Band dispersions of the 
hardcore boson in the presence of the flux as shown in Fig. 
|2]for i:> = 0.025, D = 0.1 and D = tt/S, respectively; (b), 
(d) and (f) are the corresponding bottom band curvatures for 
(a), (c) and (e). 



B. Fermionic constructions for the bosonic 
wavef unct ions 



To study the topologically ordered chiral spin state, 
we will employ the fermionic approach to construct trial 
bosonic wavefunctions . The many-body bosonic wave- 
function can be represented as follows 



1$) 



(14) 



{Xi,--- ,XNu} 



Here the sum is over all possible boson configura- 
tions \{xi}) = \{xi,--- ,xnJ) = 4r--4^rjo) 

<l>(xi,--- ,XArJ is the symmetric wavefunction. In this 
paper, we are only concerned with translationally invari- 
ant ground states. 

The many-body bosonic wave function for a product 



4 



state (PS) has the form 



N 



(15) 



i=l 



where (l)i{xi 



(j)i{xi) to maintain translational in- 



variance. / = 1, 2, 3 denotes sites within the unit ceh and 
the vector a^, j = 1, 2, is a Bravais vector for the kagome 
lattice. So the many-body wave function is labehed by 
three complex parameters ^i, 02, and 03 corresponding 
to the spin orientation on the three sites in each unit cell. 
The mean field ground state of this type is obtained by 
minimizing the average energy by varying 01, 02, and 
03. This type of spin ordered states without topological 
order is the main competing state for the ground state of 
our model. 

To construct the bosonic ground state with topological 
order, we split the the hardcore boson into two species of 
fermions 



(16) 



where ai and I3i are fermion operators which satisfy the 
constraint on every site: n^Q, = Uifs = n^^. The configu- 
ration becomes 

|{x,}) = at^/3t^---4„^/3t^JO) (17) 

and the symmetric wave function factorizes as follows 

= ^a{{x,mp{{x,)) (18) 

where ^(^{{xi}) and ^ f3{{xi}) are the antisymmetric 
fermionic wavefunctions for ai and pi. Using this frac- 
tionalization we construct two ansatz wavefunctions: the 
bosonic quantum Hall state (QHS) and the spin Hall 
state (SHS). 

The fermionic wavefunctions ^ a{{xi}) and ^/^({x^}) 
can be constructed from the mean field tight binding 
Hamiltonian 



Ha = -^eff ^ exp(z74g ) + h.c. 

{ij) 

Hp = -teff^(/3l/3,exp(z4)+h.c.) (19) 

(ij) 

The filling factors for the fermions ai^Pi are / = 1/6 per 
site, namely half per unit cell. For the QHS and SHS, we 
need the filling factor corresponding to one particle per 
unit cell which can be realized by inserting half a fiux 
quantum (0 = tt) in the original unit cell to double the 
unit cell: 



(20) 



where (j)^ and 02 are fiuxes in the kagome unit cell, see 
Fig. |2j In the presence of these fiuxes, we specify a gauge 



for Af- in this tight-binding model (19) to obtain single 



particle wavefunctions in the bottoin band: ip^Cki^Xj 



{uj = a,/3), where k is the Bloch momentum vector for 
the doubled unit cell. Thus the fermionic wavefunctions 
are the determinants of these single particle wave func- 
tions: 



"^ojiixi}) = det[V^,^(k^,: 



(21) 



where uj = /3 and i,j 
state, we set 0^ — 
we set 6? = —6^ and — 



1,2,- ,7V5. For the QHS 



^ and 02 



A/5 ■ 



for the SHS state. 



1 — — ciina — ~v^2- For the QHS, each 
fermion has the same Chern number C^j = l'^, for the 
SHS, Ca = 1 and Cf3 = -1. 

We now consider the effective theory for the QHS and 
SHS. For the QHS and SHS, fermionic excitations are 
gapped out. There is a gauge freedom for the frac- 
tionalization in Eq. (16): the gauge transformation 
ai aie'^^\ f3i /3^e~^ does not change the bosonic 
operator bi. The Hamiltonian with the gauge fluctuations 
is given by 

Ha = -^eff ^ (c^jaj exp{iAij + iaij) + h.c.^ 

(ij) 

H(3 = -teff^(/3|/3^exp(-ai^-m^^)+h.c.)(22) 

{ij) 

The non-zero Chern number for each fermion species im- 
plies that the low energy effective action for the gauge 
fields is given by 



(23) 



where ... represents higher order terms. For the QHS, 
Ca = Cf3 = l^o we obtain the low-energy effective theory 



1 



27r 



(^fiuxa^,dj,ax + (24) 



This describes the = 1/2 FQH state for bosons corre- 
sponding to the chiral spin state first introduced in Ref. 
1511 We note that although the a and P fermions have 
the same Chern number, the sign of the coupling of each 
fermion to the U{1) gauge field is opposite. Thus a 
27T fiux of creates an a fermion and annihilates a (3 
fermion. 

For the SHS, = 1, = — 1 and we obtain the 
low-energy effective theory 



>CsHS 



(25) 



Here the a and /3 fermions in the SHS have opposite 
Chern number and the sign of the coupling of the two 
fermions to the U{1) gauge field is also opposite. Thus 
a 27T fiux of creates an a fermion and a P fermion which 
corresponds to a 6 boson, i.e. it describes a spin fiip. So 
the magnetic field of the U{1) gauge field corresponds 
to the spin density. Since spin is conserved, the 
U{1) instanton is forbidden. Thus the 2+ ID U{1) gauge 
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FIG. 5: (color online) A scheme to tune the filling num- 
ber of hardcore bosons bi: the kagome lattice couples to a 
ferromagnetic substrate by the exchange interaction i^int = 
Jex • S^; we can tune the substrate magnetization (S^^) 

by an applied magnetic field. 



FIG. 4: (color online) Energy per site for D = 0.025, D = 0.1 
and D — n/S respectively. 



theory above is not confined and the U{1) gauge field 
remains gapless. As this gapless U{1) gauge field cor- 
responds to the spin density, the gapless spin den- 
sity fluctuations imply that e^^^"" spin rotation is spon- 
taneously broken. Thus the SHS is a spin XY ordered 
state. 



C. Numerical results 

For the three different states (PS, QHS and SHS), we 
can evaluate the expected energies for the bosonic model 
pi: 



mm) 



\{{x^m? 



(26) 



where we define the local energy eL{{xi}) = • 
We evaluate the energy in ( |26| ) by appropriately averag- 
ing the local energy eL{{xi\)oYeT a set of configurations 
\{xi}) distributed according to the square of the wave 
function generated with a standard varia- 

tional Monte Carlo method. 

Then we use a minimization function to optimize the 
expectation values on an 8 x 8 lattice. In Fig. |4j we plot 
the energy per site of the three states for D = 0.025, 
D = 0.1 and D = 7r/8. For D = 0.025, the PS and SHS 
have energies close to each other. The QHS has a better 
energy. As the spin-orbit coupling is increased, the SHS 
becomes worse in energy. Both the PS and QHS gain in 
energy and the PS gains much more. When D = 7r/8, 
the PS gives results close to the QHS. With small spin- 
orbit coupling {D = 0.025 and D = 0.1), the bottom 
band of the hopping Hamiltonian in ( p^O] ) is fiat and has 
a smooth curvature over the Brillouin zone. The classic 
PS cannot gain much energy through condensation of 
the lowest states. When D increases, the bottom band 
becomes more convex and the PS will gain a lot of energy 
through condensation. 

Our numerical results indicate that the topologically 
ordered QHS (the chiral spin state) is a serious candidate 
for a kagome spin system with spin-orbit coupling and 



spin polarization. This may be a realistic route for the 
discovery of new topologically ordered states in quantum 
spin systems. 



D. Practical realization 

Here we explore the possibility of obtaining a po- 
larized state {S^ ^ 0) experimentally. This can be 
achieved by applying a magnetic field, which adds a term 
Hh = —B^-Ui to the Hamiltonian (10). However, the 
exchange energy J ~ lOOmeV is usually very large, so 
experimentally accessible magnetic fields cannot polarize 
the spin to = 1/3. Hence we should find other ways 
of obtaining a large effective magnetic field. 

One way is to place the kagome lattice on a ferromag- 
netic substrate, see Fig. |5] The exchange interaction is 
^int = ^ex Xli ' ^^^^ ^ex IS the exchange coupling 
between spins on the substrate and spins Si on the 
kagome plane. The exchange coupling Jgff can be very 
large when the kagome plane structure matches that of 
the substrate perfectly. Such a large effective magnetic 
field can polarize the spin on the kagome lattice. 

A third way is to insert ferromagnetic atoms in the 
kagome system. If these ferromagnetic atoms form a fer- 
romagnetic state, the exchange interaction can also in- 
duce spin polarization on the kagome lattice. 



IV. SUMMARY 

In this paper, we study quantum spin systems on the 
kagome lattice with spin-orbit coupling and and spin 
polarization. We argue that such a system can be in 
a topologically ordered chiral spin state, a FQH state 
for bosonic spin degrees of freedom. The energy scale 
of the bosonic FQH state is of the same order as the 
spin-orbit coupling and ferromagnetism — overall much 
higher than the energy scale of FQH states in semicon- 
ductors. This result suggests exploration of topologically 
ordered states in quantum spin systems with a proper 
combination of geometric frustration, spin-orbital cou- 
pling and ferromagnetism. 

This research is supported by NSF Grant No. DMR- 
1005541 and NSFC 11074140. 
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Appendix A: Discussion of spin-orbit coupling 

The Rashba spin-orbit coupling is weak, around \D\ = 
0.025 for Herbertsimthsite Zn Cu3(OH)gCl2. This smah 
value prompts us to find other mechanisms to increase 
the strength of the spin-orbit coupling. 

In the literature, the spin-orbit coupling is also dis- 
cussed on the atomic level and with a large strength of 
coupling, e.g. around 0.2 eV and 0.4 eV for the Ad and 
5d electrons, respectively. We hope to relate the atomic 
spin-orbit coupling to the form represented in Eq. 
This has been achieved for the electron in Sr2lrOi 
and we discuss the general case below. 

A S = 1/2 electron can be found in d^, d^ and d^ 
orbtitals in transition metal cations, e.g. in Mo^^, Ir^^ 
and Cu^^ respectively. Due to the crystal field , the 
fivefold degenerate d state is split into a doublet Cg and 
a triplet t2g. d^ and d^ with a ligand octahedron and d^ 
with a ligand tetrahedron belong to t2g . The triplet t2g 
has strong spin-orbit coupling 



Hi = AL 



(Al) 



Here is the spin operator and / = 1 is the effec- 
tive angular momentum with |/f = 0) = \XY)i and 
|/f = ±1) = -^{i\XZ), ± \YZ),) (X, Y and Z are 
local axes supporting by the local octahedron or tetrahe- 
dron , e.g. see Fig. [ijb)). The strong spin-orbit coupling 
( Al ) splits t2g into two groups with effective angular mo- 
mentum J^f^ = 1/2 and J^^ = 3/2, respectively. The 
Jgff = 1/2 singlet contains a Kramers doublet: = 

^1 - i,t), + ^|o,;), and i;), = ^1 + 1,;)^ + ^|o,t).. 

Here we are concerned with the Kramers doublet labelled 
by the pseudospin = a^/2. 

There is a p bond between the cation and mediating 
oxygen atom. When the two cations and oxygen lie on a 
straight line along the bond, the values of the overlap for 
{Iz = 0\p)o and {Iz = ±l|p)o are the same. In Herbert- 
smithite, they are not along a straight line and form a 
triangle, see Fig. [ijc). As a result, different orbits have 
different overlaps with the p orbital in the oxygen {d-\p)o 



when the electron hops on the bond ri2 from ri to 

|(;, = ±l|p)o| > K?. = 0|p)o| (A2) 

For the sake of simplicity, we neglect the overlap \{lz = 
0|p)o| here and are only concerned with \{lz = ±l|p)o|- 
Then we can write the spin-orbit in terms of the pseu- 
dospin in this manner 



^ Al(rO • §(r,) 



(A3) 



Here the orbital angular momentum l(r^) can be regarded 
as the effective magnetic field on the pseudospin s(r^). It 
is very interesting that the orientation of the magnetic 
field \{Ti) varies from each site r^. When the hopping 
encloses a loop, e.g. A123 in Fig. [l] (b), an electron 
obtains a non-zero Berry phase (/)Berry related to the spin- 
orbit coupling vector D in Eq. ([T]) as follows 



^Berry = f C^r • D(r) 



(A4) 



where ft is the closed loop. 

To rotate the Kramers doublet from the local axes (X, 
Y and Z) to the global axes (x, y and z) we use 



\d;{x}) = {e^^-"^'^^U\a;{X}) 



(A5) 



where {x} 



'{X}. The hopping process on the 



bond ri2 is given as 

tl2 = ^ (taa'{^l,^2)c[^C2a' + h.C.^ (A6) 



with the hopping parameter 



^aa'(ri,r2) (e" 



l(±1|0^)(0^|±1)2 



-t{e 



-ia-ni2^ 



(A7) 



This is exactly the same as in Eq. ([T]) when D12 = ni2^. 
The Berry phase can be very large on the Kagome lattice 
resulting a very strong spin-orbit coupling effect. 
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